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Abstract. We consider two different toy cosmological models based on two fields 
(one normal scalar and one phantom) realizing the same evolution of the Bang- 
to-Rip type. One of the fields (pseudoscalar) interacts with the magnetic field 
breaking the conformal invariance of the latter. The effects of the amplification 
of cosmic magnetic fields are studied and it is shown that the presence of such 
effects can discriminate between different cosmological models realizing the same 
global evolution of the universe. 
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1. Introduction 

The discovery of the cosmic acceleration [1] and the search for dark energy responsible 
for its origin [2] have stimulated the study of different field models driving the 
cosmological evolution. Such a study usually is called the potential reconstruction 
[3], because the most typical examples of these models are those with a scalar field, 
whose potential should be found to provide a given dynamics of the universe. In the 
flat Friedmann models with a single scalar field, the form of the potential and the time 
dependence of the scalar field are uniquely determined by the evolution of the Hubble 
variable (up to a shift of the scalar field). During last years the models with two 
scalar fields have also become very popular. This is connected with the fact that some 
observations point out that the relation between the pressure and the energy density 
could be less than -1 [4]. Such equation of state arises if the matter is represented by a 
scalar field with a negative kinetic term. This field is called "phantom" [5]. Moreover, 
according to some observations [6] the universe undergoes a transition between normal 
and phantom phase. Such an effect is dubbed "phantom divide line crossing" [7]. In 
principle, the effect of phantom divide line crossing can be explained in the model with 
the only scalar field provided a special form of the potential and initial conditions is 
chosen [8] or in models with a non- minimally coupled scalar field [9]. However, the 
models with two scalar fields, one standard and one phantom, look more "natural" for 
the description of the phantom divide line crossing [10, 11, 12]$. 

In our preceding paper [12] we have studied the procedure of reconstruction of the 
potential in two-field models. It was shown that there exists a huge variety of potentials 
and time dependences of the fields realizing the same cosmological evolution. Some 
concrete examples were considered, corresponding to the evolution beginning with the 
standard Big Bang singularity and ending in the Big Rip singularity [13]. 

One can ask oneself: what is the sense of studying different potentials and scalar 
field dynamics if they imply the same cosmological evolution? The point is that the 
scalar and phantom field can interact with other fields and influence not only the 
global cosmological evolution but also other observable quantities. 

One of the possible effects of the presence of normal and phantom fields could 
be their influence on the dynamics of cosmic magnetic fields. The problem of the 
origin and of possible amplification of cosmic magnetic fields is widely discussed 
in the literature [15]. In particular, the origin of such fields can be attributed to 
primordial quantum fluctuations [16] and their further evolution can be influenced by 
hypothetic interaction with pseudoscalar fields breaking the conformal invariance of 
the electromagnetic field [17, 18]. In the present paper we consider the evolution 
of magnetic fields created as a result of quantum fluctuations, undergoing the 
inflationary period with unbroken conformal invariance and beginning the interaction 
with pseudoscalar or pseudophantom fields after exiting the inflation and entering the 
Big Bang expansion stage, which is a part of the Bang-to-Rip scenario described in 
the preceding paper [12]. We shall use different field realizations of this scenario and 
shall see how the dynamics of the field with negative parity influences the dynamics 
of cosmic magnetic fields. To our knowledge the possible influence of the two-field 
dynamics, (when one of two (pseudo) scalar fields is a phantom one) on the magnetic 
fields was not yet discussed in the literature. 

Speaking of cosmic magentic fields we mean the large-scale galactic, intergalactic 

| These two fields may have their origin in spontaneous breaking of primordial symmetries of moduli 
space with complex potential [11]. 
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or super-cluster magnetic fields of order from 10 -6 G to 10~ n G with correlation 
from 100 kpc to several Mpc to the extent that they are originated from scalar and, 
possibly gauge field fluctuations after exiting the inflation. Their seeds may well have 
10~ 18 - 10- 27 G or less (see [15]). 

The structure of the paper is as follows: in Sec. 2 we recall the Bang-to-Rip 
scenario and describe some examples of different dynamics of scalar and phantom 
fields; in Sec. 3 we introduce the interaction of the fields (phantom or normal) with an 
electromagnetic field and write down the corresponding equations of motion; in Sec. 4 
we describe the numerical simulations of the evolution of magnetic fields and present 
the results of these simulations; Sec. 5 is devoted to concluding remarks. 



2. Cosmological evolution and (pseudo)-scalar fields 

We shall consider a spatially flat Friedmann universe with the metric 

ds 2 = dt 2 -a 2 (t)dl 2 . (1) 

Here the spatial distance element dl refers to the so called comoving coordinates. The 
physical distance is obtained by multiplying dl by the cosmological radius a(t). We 
would like to consider the cosmological evolution characterized by the following time 
dependence of the Hubble variable h(t) = a/a, where "dot" denotes the differentiation 
with respect to the cosmic time t: 

This scenario could be called "Bang-to-Rip" [12] because it is obvious that at small 
values of t the universe expands according to power law: a <~ i 1 / 3 while at t — > tn 
the Hubble variable explodes and one encounters the typical Big Rip type singularity. 
(The factor one third in (2) was chosen for calculation simplicity). In our preceding 
paper [12]) we considered a class of two-field cosmological models, where one field was 
a standard scalar field <j>, while the other was a phantom one £. The energy density 
of the system of these two interacting fields is 

e=\j> 2 -\i 2 + V{<f>,d). (3) 

Analyzing the Friedmann equation § 

h 2 = e, (4) 

we have shown that in contrast to models with one scalar field, there is huge variety 
of potentials V(<p, £) realizing a given evolution. Moreover, besides the freedom in the 
choice of the potential, one can choose different dynamics of the fields 4>{t) and £(t) 
realizing the given evolution. We have studied in [12] some particular exactly solvable 
examples of forms of potentials and time dependences of fields <j> and £ providing the 
evolution (2). Here we shall present and apply some of them. Consider the potential || 

^ = ^2 cosh 6 (-30/4) exp(3V20- (5) 

R 

§ We use the following system of units K = 1, c = 1 and 8nG = 3. In this system the Planck mass 
mp, the Planck length Ip and the Planck time ip are equal to 1. Then when we need to make the 
transition to the "normal", say, cgs units, we should simply express the Planck units in terms of the 
cgs units. In all that follows we tacitly assume that all our units are normalized by the proper Planck 
units. Thus, the scalar field entering as an argument into the dimensionless expressions should be 
divided by the factor yjmp /t p . 

|| The expression for the potential should be multiplied by the factor mp /(dp). 
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and the fields 



4>{t) = -^arctanhy-^— -, (6) 



(Here the expressions for the fields <p(t) and should be multiplied by \Jmp/tp. 
For the relation between Planck units and cgs ones see e.g. [19] If we would like to 
substitute one of these two fields by the pseudoscalar field, conserving the correct 
parity of the potential, we can choose only the field <j> because the potential Vj is even 
with respect to <fi, but not with respect to £. In what follows we shall call the model 
with the potential (5), the pseudoscalar field (6) and the scalar phantom (7) "model J". 
Consider another potential 

M£> 4>) = sinh 2 (3£/4) cosh 2 (3£/4) exp(-3\/2». (8) 

with the fields 

= ^arctanh^J. (10) 

This potential is even with respect to the field £. Hence our model II is based on the 
potential (8), the pseudophantom field (10) and the scalar field (9). They will be the 
fields with the negative parity which couple to the magnetic field. 

3. Post-inflationary evolution of a magnetic field interacting with a 
pseudo-scalar or pseudo-phantom fields 

The action of an electromagnetic field interacting with a pseudoscalar or 
pseudophantom field <f> is 

S = -\J d 4 x^g-(F^F^ + a<f>F^F^), (11) 

where a is an interaction constant and the dual electromagnetic tensor F^ v is defined 

as 

pnu = Ie^F^, (12) 

where 

J? = ^rT7, f wvvpo = \ f P-vpa fio-i 



with the standard Levi-Civita symbol 

£\ivpa = £[nvpo-] i £0123 = +1- (14) 

Variating the action (11) with respect to the field we obtain the field equations 

V„F^ = -ad^F^, (15) 

V^F"" = 0. (16) 
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The Klein-Gordon equation for the pseudoscalar field is 

dV 



d(j> 



(17) 



The Klein-Gordon equation for the pseudophantom field (which is the one that couples 
with the magnetic field in the model 77) differs from equation (17) by change of sign 
in front of the kinetic term. In what follows we shall neglect the influence of magnetic 
fields on the cosmological evolution, i.e. we will discard the electromagnetic coupling 
in equation (17). 

If one wants to rewrite these formulae in terms of the three-dimensional quantities 
(i.e. the electric and magnetic fields) one can find the expression of the electromagnetic 
tensor in a generic curved background, starting from a locally flat reference frame — 
in which it is well known the relation between electromagnetic fields and F — and 
using a coordinate transformation. It is easy to see that we have, for the metric (1): 
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The field equations (15), (16) and (17) rewritten in terms of E and B become 

V • E= -ctV(f) ■ B. 

d (a 2 E) - V x (aB) = -a[d c/>(a 2 B) - W<f> x (aE)}, 
d (a 2 B) - V x (aE) = 0, 

V • B = 0. 



(19a) 
(196) 
(19c) 
(19d) 



For a spatially homogeneous pseudoscalar field equations (19a) and (196) look 

like 

V • E = (20a) 

d (a 2 E) - V x (aB) = —ad 4>{a 2 B). (206) 

Taking the curl of (206) and substituting into it the value of E from (19c) we obtain 

A^(a 2 B) a. 



d£{a 2 B) + h(t)d {a 2 B) 



-do^V x (a 2 B) = 0, (21) 



a' a 

where A^ 3 ) stands for the three-dimensional Euclidean Laplacian operator. 
Let us introduce 

F(x, t) = a 2 {t)B(x, t) 

and its Fourier transform 

1 



F(k,t) = 



(2tt) 3 / 2 



e- tk - x F(x,t)d 6 x. 



(22) 
(23) 



Here the field B is an observable magnetic field entering into the expression for the 
Lorentz force. The field equation for F(k, t) is 

2 



F(k,t) + h(t)F(k,t) + 



— (j>kx 
a 



F(k,t) = 0, 



(24) 
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where "dot" means the time derivative. This last equation can be further simplified: 
assuming k = (k, 0, 0) and denning the functions F± = (F 2 ± iFs)/y/2 one arrives to 



F± + hF± 



a J a 



F± = 0, (25) 



where we have omitted the arguments k and t. 

Assuming that the electromagnetic field has a quantum origin (as all the fields 
in the cosmology of the early universe [14]) the modes of this field are represented 
by harmonic oscillators. Considering their vacuum fluctuations responsible for their 
birth we can neglect the small breakdown of the conformal symmetry and treat them 
as free. In conformal coordinates (77, x) such that the Friedmann metric has the form 

ds 2 = a 2 {r]){dr] 2 - dl 2 ), (26) 

the electromagnetic potential Ai with the gauge choice Aq — 0, djA 3 = satisfies the 
standard harmonic oscillator equation of motion 

A\ + k 2 A t = 0. (27) 



Hence the initial amplitude of the field Ai behaves as Ai — 1/V2fc, while the initial 
amplitude of the functions F is y/k/2. The evolution of the field F during the 
inflationary period was described in [18], where it was shown that the growing solution 
at the end of inflation is amplified by some factor depending on the intensity of the 
interaction between the pseudoscalar field and magnetic field. 

Here we are interested in the evolution of the magnetic field interacting with the 
pseudoscalar field after inflation, where our hypothetic Bang-to-Rip scenario takes 
place. More precisely, we would like to see how different types of scalar-pseudoscalar 
potentials and field dynamics providing the same cosmological evolution could be 
distinguished by their influence on the evolution of the magnetic field. We do not take 
into account the breaking of the conformal invariance during the inflationary stage 
and all the effects connected with this breakdown will be revealed only after the end 
of inflation and the beginning of the Bang-to-Rip evolution. This beginning is such 
that the value of the Hubble parameter, characterizing this evolution is equal to that 
of the inflation, i.e. 

M*<>) = q , ,l R - ; = inflation * lO^S" 1 . (28) 
oto{lR — to) 

In turn, this implies that we begin evolution at the time moment of the order of 10 -33 s. 

We shall consider both the components F+ and F- and we shall dwell on the 
scenarios I and II described in the preceding section. Anyway, our assumption 
regarding the initial conditions for equation (25) can be easily modified in order 
to account for the previous possible amplification of primordial magnetic fields as 
was discussed by [18]. Thus, all estimates for the numerical values of the magnetic 
fields in today's universe should be multiplied by some factor corresponding to the 
amplification of the magnetic field during the inflationary stage. Hence, our results 
refer more to differences between various models of a post inflationary evolution than of 
the real present values of magnetic fields, whose amplification might be also combined 
effect of different mechanism [15, 17, 16, 18]. 
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4. Generation of magnetic fields: numerical results 

In this section we present the results of numerical simulations, for the two models I 
and II introduced in Sec. 2. In our models, the equation of motion for the modes 
F±(k,t) (25) readst 



F ± + -. tRF± 



3t(t R - 1) 



F± = 0,(29) 



where <& stands for the scalar field <f> in the model I and for the phantom £ in the 
model II, so that 



$I ^ = 3^' ^^^svw^ry (30) 

Equation (29) is solved for different values of the wave number k and the coupling 
parameter a. (The parameter a has the dimensionality inverse with respect to that 
of the scalar field; the wave number k has the dimensionality of inverse length; the 
time tn = 10 17 s). Qualitatively we remark that in (25) the coupling term influence 
becomes negligible after some critical period. After that the magnetic fields in our 
different scenarios evolve as if the parameter a in (25) had been put equal to zero. 
Indeed, it can be easily seen that the interaction term vanishes with the growth of the 
cosmological radius a. Then the distinction between the two models is to be searched 
in the early time behavior of the field evolution. 

Noting that in both our models I and II the time derivative 4> is positive, by 
inspection of the linear term in equation (25) we expect the amplification to be mainly 
given for the mode F_ provided the positive sign for a is chosen; so we will restrict our 
attention on F_ . We can also argue that the relative strength of the last two terms in 
the left-hand side of (29) is crucial for determining the behavior of the solution: when 
the coupling term prevails (we remark that we are talking about F_ so this term is 
negative in our models) then we expect an amplification, while when the first term 
dominates we expect an oscillatory behavior. For future reference it is convenient to 
define 

A(t;*)^-^ = 4(^y /3 , (3D 



a(t)a§ a$ V t 

which is just the ratio between the last two terms in the left-hand side of equation (29). 

Indeed our numerical simulations confirm these predictions. Let us consider the 
model I with a = l(tp /mp) 1 / 2 and k = lO -55 /^ 1 , where lp is the Planck length. Such 
a value of the wave number k corresponds to the wave length of lkpc at the present 
moment. We obtain an early-time amplification of about 2 orders of magnitude, with 
the subsequent oscillatory decay. Notice that the parameter A in this model at the 
beginning is very small: this corresponds to the dominance of the term proportional 
to $ and, hence, to the amplification of the field F_. At the time scale of the order 
of 10 51 ip, where tp is the Planck time, this regime turns to that with big values of A 
where the influence of the term proportional to $ is negligible. 

For the same choice of the parameters a and k in the model 77 the amplification 
is absent. 

K The reader can easily verify that this equation is obtained imposing the normalization a| to day = 1, 
where today-time is taken to be near the crossing of the phantom divide line, i.e. at t ~ tft/2. 
This implies in turn that at the beginning of the "Bang-to-Rip" evolution the cosmological radius is 
a(t ) ~ 1CT 17 . 
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Figure 1. Plot of the ratio A for the model /, with the parameter choice 
a = l(tp /mp) 1 / 2 , k = lO -55 ^ 1 . It can be easily seen that at a time scale 
of order 10 51 tp the ratio becomes greater than 1. 




1.x 10^ 2.x 10^ 3.x 10^ 4.X1CT 

Figure 2. Plot of the time evolution of the absolute value of the complex field F 
(given in Planck units) in model / with the parameter choice a = lftp/rap) 1 / 2 , 
k = lO -55 ^ 1 . The behavior, as said above in the text, consists in an amplification 
till a time of order 10 51 tp, after which the oscillations begin. 
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2.x 1(T 4.x 1(T 6.x KT 8.xl(T 1.x 10" 

Figure 3. Plot of the time evolution of the absolute value of the complex field F 
(given in Planck units) in model / with the parameter choice a = 100(tp/mp) 1 / 2 , 
k = 10~ 55 lp . The behavior, consists in an amplification till a time of order 
10 54 tp, after which the oscillations begin. 



In Figure 1 we present the time dependence of the function A for the model I 
for the values of a = l^p/mp) 1 / 2 and k — lO -55 ^ 1 chosen above. The Figure 2 
manifests the amplification of the magnetic field in the model /. 

Naturally the effect of amplification of the magnetic field grows with the coupling 
constant a and diminishes when the wave number k increases. In Figure 3 we display 
the results for the case of a = 100(ip/mp) 1 / 2 , which is admittedly extreme and 
possibly non realistic, but good for illustrative purposes. Here the amplification is 
more evident and extends for a longer time period. 

Let us try to make some estimates of the cosmic magnetic fields in the universe 
today, using the correlation functions. The correlation function for the variable F is 
defined as the quantum vacuum average 

G l3 (t,x-y) = (0|^(M)F^,y)|0) (32) 

and can be rewritten as 

/d 3 k 

K e -iHz-yl F .(t >k)F *(t >k) . (33) 

Integrating over the angles, we come to 

G^-y) = J ^ ^S^ mmm. (34) 

To estimate the integral (34) we notice that the main contribution to it comes from 
the region where k « l/\x — y\ (see e.g. [19]) and it is of order 

^(M/£)l 2 , (35) 

where L = \x — y\. In this estimation the amplification factor is 
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where the subscript i is not present since we have taken the trace over polarizations. 

Now we are in a position to give numerical values for the magnetic fields at 
different scales in the model / for different values of the coupling parameter a. These 
values (see Table 1) correspond to three values of the coupling parameter a (1,10 and 
100 (tp/mp) 1 / 2 ) + and to two spatial scales L determined by the values of the wave 
number k. We do not impose some physical restrictions on the value of a. It is easy 
to see that the increase of a implies the growth of the value of the magnetic field B. 
We also shall consider a rather large value of the wave number k corresponding to 
the physical wavelength ~ 10~ 2 pc at the present moment when a = 1. While this 
scale looks too small for the description of large-scale cosmic magentic fields, we use 
it for illustrative purposes to underline the strong dependence of the amplification of 
magnetic field on the corresponing wavelength values. 

Let us stress once again that we ignored the effects of amplification of the magnetic 
fields during inflation to focus on seizable effects during evolution. 



Table 1. The Table displays the values of the magnetic field B 
corresponding to the chosen values of a and k. The length L 
refers to the present moment when a = 1. 

















a = lyjtp/mp 


a. = 10yjtp/mp 


a = I00y/tp/mp 


k = lO- 55 ^ 1 (L = 


lkpc) 


B ~ 10- 67 G 


B ~ 10- 60 G 


B ~ 100G 


k = lO" 50 ^ 1 (L = 


10"V) 


B ~ 10~ 55 G 


B ~ 10~ 49 G 


B ~ 10 13 G 



Finally notice that our quantum "initial" conditions correspond to physical 
magnetic fields which for presented values of k are Bi n ~ k 2 /a 2 is equal to 10~ 34 G 
for k = lO- 55 /^ 1 and B in - 10- 24 G for k = lO" 50 /^ 1 . 

5. Conclusion 

We have seen that the evolution of the cosmic magnetic fields interacting with a 
pseudoscalar (pseudophantom) field is quite sensitive to the concrete form of the 
dynamics of this field in two-field models where different scalar field dynamics and 
potentials realize the same cosmological evolution. 

We confirm the sensitivity of the evolution of the magnetic field with respect to 
its helicity given the sign of the coupling constant a and that the $ is a monotonic 
function of time (as it is really so in our models). We give also some numerical 
estimates of the actual magnetic fields up to the factor of amplification of such fields 
during the inflationary period. The toy model of the Bang-to- Rip evolution studied in 
this paper, cannot be regarded as the only responsible for the amplification of cosmic 
magnetic fields implying their present observable values. It rather complements some 
other mechanisms acting before. However, the difference between cosmic magnetic 

+ It is useful to remark that at this length scales values of a less than 1 make the effect of the coupling 
of the magnetic field with the pseudoscalar field negligible. 
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fields arising in various models (giving the same expansion law after the inflation) is 
essential. It may provide a discriminating test for such models. 

Naturally, the influence of the interaction between a pseudoscalar (phantom) field 
and a cosmic magnetic field on the dynamics of the latter depends on the velocity of 
change with time of the former. The larger is the time derivative of the pseudoscalar 
field, the more intensive is the growth of the magnetic field (remember that the 
evolution of the scalar field is monotonic). The results of numerical calculations 
illustrated in Sec. 4 confirm these qualitative considerations. Moreover, one can 
see that there exists a certain range of values of the wave number k (and hence of 
the corresponding wavelengths of the magnetic field) where the effect is stronger. 
Indeed, if k is too small the interaction term is small as well and the evolution of the 
magnetic field is damped. On the other hand if the wave number k is too large the 
interaction term is small compared to the oscillatory term, proportional as usual to k 2 
and the evolution has practically oscillatory character. Thus, the study of interplay 
between the dynamics of global scalar fields providing the cosmological evolution and 
the magnetic fields looks promising. This interplay has another interesting aspect. The 
pseudoscalar-electromagnetic field interaction can imply a conversion of the photons 
into axions. Such an effect can cause the observable dimming of supenovae. While 
it was shown [20] that this effect cannot mimic the cosmic acceleration, it could 
nevertheless mimic the dark energy fluid with a phantom equation of state [21]. Thus, 
the interrelation amongst an electromagnetic field, a scalar and a phantom field can 
reveal some surprises. 
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